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We study analytically a model where particles with a hard-core repulsion diffuse on 
a finite one-dimensional lattice with space-dependent, asymmetric hopping rates. 
The system dynamics are given by the U g [SU(2)]-symmetric Hamiltonian of a 
generalized anisotropic Heisenberg antiferromagnet. Exploiting this symmetry we 
derive exact expressions for various correlation functions. We discuss the density 
profile and the two-point function and compute the correlation length £ s as well as 
the correlation time £j . The dynamics of the density and the correlations are shown 
to be governed by the energy gaps of a one-particle system. For large systems £ s 
and £t depend only on the asymmetry. For small asymmetry one finds £j ~ £;? 
indicating a dynamical exponent z = 2 as for symmetric diffusion. 
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Recently there has been much interest in stochastic processes the time 
evolution of which can be mapped to Hamiltonians of quantum lattice mod- 
els. A particularly interesting class of systems (comprising many interesting 
equilibrium and non-equilibrium models) are those where the dynamics are 
mapped to integrable quantum chains in one dimension [[L| - || or to SU(2)- 
symmetric quantum Hamiltonians in any dimension ||. In these cases the 
symmetries of the models can be used to determine critical exponents, cor- 
relation functions and other quantities of interest. Among the models which 
can be studied in this way are diffusive lattice gases (which are related to 
growth models), diffusion-limited chemical reactions and vertex models. 

In this letter we study a driven lattice gas of particles on a one- 
dimensional chain of length L. We assume a hard core repulsion preventing 
the occupation of each site by more than one particle. Particles hop stochasti- 
cally to their right and left nearest neighbours with space-dependent hopping 
rates. We consider reflective boundary conditions, i.e., particles cannot enter 
or leave the system. We will show that if the driving force (represented by 
an asymmetry in left- and right-hopping probabilities) is spatially constant 
the system has a U g [SU(2)] symmetry even if the particle mobility has an 
arbitrary space dependence. We will use this symmetry to derive exact re- 
sults for the correlation functions even for finite systems. One should keep 
in mind that even though we are studying a specific model, our method and 
some of the results can be generalized to other systems. This aspect will be 
briefly discussed at the end of this letter. 

In what follows we will study one- dimensional systems described by oc- 
cupation numbers n = {rij} where j is the number of a lattice site und runs 
from 1 to L. Their dynamics are given by a master equation for the proba- 
bility distribution F(n, t) which can be mapped onto a Schrodinger equation 
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H-0: 

d t \F(n,t)) = H\F(n,t)) . (1) 

We will focus on a system of particles with hard core repulsion diffusing 
on a lattice with position dependent, asymmetric hopping rates. It shall 
be subjected to reflecting (sometimes called free) boundary conditions. The 
Hamiltonian is then given by: 

L-l 

3=1 
L-l 

= Y,\[*XiPi u i ■ ( 3 ) 

3=1 

Here c] creates a particle at site j and Cj annihilates a particle at site j. 
They obey Pauli commutation relations, aij is the probability rate for a 
single particle jump from site j + 1 to site j and (3j is the rate for a jump 
from j to j + 1. The logarithm / 3 - = In ^- of the ratio of the hopping rates 
corresponds to a " driving force" while the product aj/3j is a measure for the 
mobility of the particles. Here we assume the force to be independent of the 
position, i.e., 

e ^ = ~q~ = ^ ^ 3 with q real . (4) 

Pj 

q = 1 corresponds to diffusion without a driving force (symmetric diffusion). 

We would like to point out that the model is related to an integrable quan- 
tum chain. Introducing operators {aj, crj, crj} = {?~ J c] + q^Cj, — i(g _J c] — 
5 3 'cj), 1 — 2cjCj}, which commute like Pauli- matrices, the Hamiltonian reads: 



(5) 

This is the Hamiltonian of a generalized Heisenberg antiferromagnet. It is 
the sum of generators Uj of a Temperley-Lieb algebra and it is symmetric 



under the action of the quantum group U 9 [SU(2)] || [1Q| . 

The generators of the quantum group U g [SU(2)] are given in terms of 
creation and annihilation operators as: 

S+ = E h{q) , S~ = jr b){q) , S* = + \) (6) 



j=l j=l j=l 



with 



b)( q ) = gEtlK-^ctg-Et^!^- 1 ) , ( 7 ) 
b^q) = g EU 1 i»*c ig -ELi+i n * (8) 



and rij = CjCj. They satisfy the following relations: 



= [2S z ] q and [5",5 ± ]_ = ±5 ± (9) 



g -Q 



where [x], 

It is easy to check that each term Uj in (||) commutes with and S* 2 . 
Hence 

[H, S' ± ]_ = (10) 
= . (11) 

Since S z = L/2 - EL 1 rij eq. ( |TTD states that the particle number is con- 



served, which is obvious for a purely diffusive system. Eq. (|10l), however, 
has the nontrivial consequence that we can construct all the zero energy 
eigenstates out of the zero particle state | ). One obtains 



H\N) = with | N) = J^(S-) N \0) (12) 
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and (N\H = with ( JV | = ( \(S + ) N . (13) 



Here [m] q \ = [l] q [2] q ...[m] q . The left eigenstates assign equal weight to any N- 
particle configuration. Hence, averaging over an iV-particle state is performed 
by multiplying ( iV | from the left. The normalized iV-particle states 

\N) norm = ( L X X \N) = (^YV iL+1) £\-^U^\n) (14) 
are the steady states of the system and satisfy ( N \ N ) = 1. Here \J^J 



[l-n] \[N] ~. > an< ^ ^ ne u PP er index at the sum means that the summation runs 
only over states with total particle number N. 

Next we will discuss the properties of the steady state in more detail. 
A configuration with the A^-particles located at sites {k%, &2, ...kj^} has a 
probability proportional to q~ 2 ^i=i kl if the system is in the steady state 
(|TJJ). Note that these states obey a datailed balance condition. Two special 
cases are easily understood. For q = 1, i.e., for symmetric diffusion, every 
configuration has the same probability in accordance with what is known. In 
the limit q — >• we find a configuration where the particles occupy the sites 
{L — N + 1, L — N + 2, L} with probability 1, which is plausible. The case 
q — > oo is equivalent to q — > since simultaneously replacing q by q^ 1 and 
each position x by L + 1 — x leaves the system invariant. 

For the calculation of averaged quantities in a state fli~4] ) the following 
relations are useful: 

c*\N) = q L+1 - 2x (l-n x )\N-l) (15) 
(N\ct = (N\(l-n x ) . (16) 
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From this we find a recursion for one-time correlation functions: 



[N] q q L+1 ~ 2xi 

( N \n xi ...n xi \ N ) norm = - — - - ( N-l |n a . 1 ...n a . l _ 1 (l-7i a . I )| N-l ) norm 

(17) 

Of particular interest is the density profile pn{x) = ( N \n x \ N ) noTm in 
the iV-particle steady state. Solving the recursion fll7D and inserting pi(x) = 
q L+1 ~ 2x /[L] q we find the exact expression: 



-1 N-l 



pn(*)=( n j2(-^) N - k+i ^ k){L+i - 2x) : • (is) 

\ / q k=0 \ / q 

The recursion relation QT7D simplifies for large and L. The following 
approximations are then derived (We assume q < 1.): 

p N {L- N + l-r) = q 2r + o(q 4r ) w e~ r/ ^ (19) 
p N (L-N + r) = 1 - q 2r + o(q Ar ) w 1 - e- r/ ^ s for r = 1, 2, ... (20) 

with the decay length 

& = ■ ( 21 ) 

From ( ]TD| ) and (^Dj) the following symmetry becomes obvious: 

p N (L _JV + l-r)wl - pjv(L - iV + r) . (22) 

In this approximation the density profile is invariant with respect to the 
scale transformation N — > bN , L — > bL , r —> br , £ s — > 6^ s . Rescaling 
only iV, L, r and keeping g fixed we obtain a theta-function for the density 
profile in the limit b — > oo. I.e., a given system observed on a large scale 
gives p{x) ~ 9{x — L + AT). Using a higher resolution we find this function to 
be smeared. The decay length £ s is a measure for the width of the transition 
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zone seperating the high density region with p — 1 from the zero density 
domain. Its center (where p ~ 1/2) is located between sites L — N and 
L — N + 1. Note that £ s does not depend on iV or L. 

Another important quantity is the connected two-point correlation which 
is defined as 

C(x, y,t) = { n x (t)n y (0) ) N - p N {x)p N {y) . (23) 

where averaging is performed over the iV-particle steady state. 

For the one-time correlation function we find an exact expression solving 
the recursion relation (|i~7|): 

C{x, y, 0) = ^ 2y _ ^ 2x pN{x)p N {y) ■ (24) 

It is interesting to note that as in a similar diffusion model with open bound- 



ary conditions, where particles are injected and removed \\ 1| - [13], the two- 
point function can be expressed in terms of one-point functions. The two- 
point function for a large system has a finite amplitude only in the transition 
zone of the density profile. Inserting the expressions (|19|) and fl2"0|) it turns 
out that the correlation length is identical to the decay length £ s (see eq. 
(M)) of the density. 



In order to understand the dynamics of the system we study the time- 
dependent correlation function (^). Using the results for Pat(x) it is easy to 
see that like the one-time correlator it is nonzero only in the transition zone 
of the density profile. A more detailed analysis starts with the explicit form 
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following from (p3|): 



C{x,y,t) = fy (N\n x e Ht n y \N)-fy ( N \ n x \ N ) ( N \ n y \ N ) . 

(25) 

Inserting a complete set of energy eigenstates, defined by H | k ) = e k \ k ) 
and ( k | H = ( k \ e k we get: 

with A*(a:) - ( JV | n x \ k ) , B k {y) = {k\n y \N) (27) 



where the prime means that the steady state is excluded from the summa- 



tion. (Its contribution cancels with the second term in (|25).) Because of 
the U 9 [SU(2)] symmetry the energy eigenstates are eigenstates of S z and the 
Casimir operator S 2 = 2S + S~ + [S^JS^ + l] q . Hence, the energy eigen- 
states | k ) can be classified by their spin-componenet S z and by their 
total spin S. Since for real q the representations of U g [SU(2)] are isomor- 
phic to the ones of SU(2) |10[ the selection rules for matrix elements are 
the same for both symmetries. Consequently, the only nonvanishing ma- 
trix elements fl27D are those with states | k ) which have S = L/2 — 1 and 
S z = L/2 — N 0. These states are given by | N, k ) = ^(S-)^ 1 ] 1, k ) 
, ( N, k | = ( l,k \(S + ) N ~ 1 , respectively, where we have introduced 
one-particle eigenstates as: 



H | l,k ) = e k | l,ife ) and ( 1, k \ H = { 1, k \ e k . (28) 

Hence, the time dependence of the correlation functions is determined by 
the eigenenergies of the one-particle system, which are much easier to cal- 
culate than many-particle energies. The matrix elements (|27|) reduce to the 
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one-particle matrix elements A k (x) = ( 1 |(5+) iV - 1 n x (S^- 1 | 1, k)/([N} q \) 2 
and B k (y) = (l,k\(S+) N - 1 n y {S')^ 1 | 1 )/([iV] 9 !) 2 . The exact form of the 
correlation function can be obtained by evaluating them which (if the eigen- 
states 1 1, k) are known) is a tedious but straightforward task. 

The time evolution of the density profile can be studied in the same way 
as the correlation functions. If the system starts with an arbitrary initial N- 
particle state | F(0) ) the density is p(x,t) — ( N \ n x e Ht \ F(0) ). It decays 
to its stationary value ( Jl8| ) on time scales given by the one-particle energies 
(f2~8l). Expanding | F(0) ) in eigenstates of H one finds that contributions to 
the time dependence of p(x, t) come from those excitations e k for which the 
matrix elements ( 1, k \ (S + ) N ~ 1 \ F(0) ) do not vanish. 

Up to this point our discussion was completely general and applies actu- 
ally to any system of exclusive particles with an U g [SU(2)]-symmetric time 
evolution operator. In what follow we will discuss the system (||D with ho- 
mogeneous jump rates <x,- = q , [3j = q^ 1 Wj. For this case we find the 
eigenenergies: 

2n 

e k = 2cos(— k) - {q + q- 1 ) < for k = 1, 2, ...L - 1 . (29) 

Li 

The long-time behaviour of the correlation function (^) (and of the density 
p{x,t) if the initial state is not orthogonal to the state | N, k\ )) is governed 
by the lowest excitation e\. Consequently, we can identify 

6 = -q 1 = {q + q- 1 ~2 cos(^)}- 1 ^ {q + q' 1 - 2}~ l (30) 

with the correlation time of the system. Note that £t does not depend on L 
or N if L 3> 1. We would like to emphasize that for (g =^ 1) ^ remains finite 
for L — > oo, as opposed to the correlation time for asymmetric diffusion with 
periodic boundary conditions which diverges as L 3//2 ||. The significance of 
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this remark lies in the fact that it demonstrates the impact of a boundary 
term which alone is sufficient to change the gap structure of the Hamilto- 
nian. Also for symmetric diffusion £ t diverges, but as L 2 , independent of the 
boundary term. If the asymmetry is small, i.e., if q — 1 -C 1, one derives from 
the expressions (|21| ) and (|30| ) that correlation time and correlation length are 
related as £t ~ ^ indicating a dynamical exponent z — 2. 



Finally, we briefly discuss possible generalizations of our results and meth- 
ods to other Uq[SU(2)]-symmetric models. Instead of the Hamiltonian (g) 
one could, for example, study an asymmetric 6-vertex model with space- 
dependent vertex weights defined by the diagonal-to-diagonal transfer ma- 



trix T = (n^i %)(nt=i7Vi) with Tj = 1 — ^~P jUj 0, |, m. Here 
the hopping probabilities ctj and f3j parametrize the vertex weights and T 
defines a discrete time evolution with a parallel updating mechanism. In the 
mapping to the vertex model the density correlations discussed here corre- 
spond to arrow correlations in the plane. Other models of interest for which 
our methods could successfully be applied are the g-deformed spin-1 Heisen- 



berg chain jRJ describing reaction-diffusion processes involving two different 
kinds of exclusive particles or some of the models discussed in [|J. Of course 
it also possible to generalize our procedure to models with other g-deformed 
symmetries. A systematic approach to the study of correlation functions of 
stochastic systems using symmetries of the type discussed here will be pre- 
sented elsewhere. 
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